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We consider hairy black hole solutions of Einstein- Yang- Mills-Dilaton theory, coupled to a Gauss- 
Bonnet curvature term, and we study their stability under small, spacetime-dependent perturba- 
tions. We demonstrate that the stringy corrections do not remove the sphaleronic instabilities of 
the coloured black holes with the number of unstable modes being equal to the number of nodes of 
f~») ' the background gauge function. In the gravitational sector, and in the limit of an infinitely large 

f^ , horizon, the coloured black holes are also found to be unstable. Similar behaviour is exhibited by 

f^ ■ the magnetically charged black holes while the bulk of the neutral black holes are proven to be 

O^ ' stable under small, gauge-dependent perturbations. Finally, the electrically charged black holes are 

(-H . found to be characterized only by the existence of a gravitational sector of perturbations. As in the 

jrt ' case of neutral black holes, we demonstrate that for the bulk of electrically charged black holes no 

unstable modes arise in this sector. 
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0\ ■ I. INTRODUCTION 

o 

f^ ' To date a large number of hairy black hole geometries have been found, in a diverse range of matter models 

^-H \ coupled to various theories of gravity (for a recent review of some aspects, see m). Stable examples of black holes 

^^ 1 with hair are of particular interest for their physical relevance, and also in order to investigate the effects of hair on 

5^ ' quantum processes connected with black holes. However, many of the hairy black holes currently known are unstable, 

O ■ particularly those involving non-Abelian gauge fields ^^, where the instability is topological in nature and similar 

O"^ to that of the flat space sphaleron. The only exceptions to the above rule are the black hole solutions found in the 

Vh i framework of the Einstein-Skyrme theory B and the magnetically charged, non-Abelian black holes in the limit of 

f^' infinitely strong coupling of the Higgs field flG] or in the presence of a negative cosmological constant [|ll| . The limited 

^ ' number of stable black holes known so far makes the quest for new stable solutions both immediate and challenging. 

''K^ ' The superstring effective action at low energies, that follows from the compactification of the heterotic superstring 

''^ , theory, provides us with a generalized theory of gravity which is an excellent framework for the study of black holes. 

d ' The theory contains, apart from the usual Einstein term, a number of scalar fields, the dilaton, axion and moduli 

fields, coupled to higher-derivative gravitational terms as well as non-Abelian gauge fields. Here, we are particularly 

interested in the theory that describes the coupling of the dilaton field to the one-loop Gauss-Bonnet curvature term, 

and can also include a non-Abelian gauge field. These models are known to possess black hole solutions with hair, 



both with |12|-|l4| and without |15 1^ the gauge field. In the absence of the gauge field, the so-called dilatonic hairy 
black holes found in [|l5| were proven to be linearly stable ll^ with their stability being in accordance with their 
interpretation as a generalization of the Schwarzschild black hole in the framework of string theory. This family of 
stable black holes corresponds to the upper branch of the neutral black hole solutions found later in ||I3] and whose 
relative stability with respect to a second unstable branch of solutions was demonstrated by the use of catastrophe 
theory (l8|. However, the question of the stability of the corresponding coloured black holes, that arise in the presence 
of the non-Abelian gauge field, still remains open. In this article we investigate whether the stability of the dilatonic 
black holes extends to the case where the dilaton is also coupled to the gauge field. We will be particularly interested 
in stringy coloured black holes, and it will be shown that they, like their non-stringy counterparts, are topologically 
unstable. In addition, conclusions on the stability of the dilatonic black holes under small, gauge perturbations as 
well as that of the magnetically and electrically charged black holes will also be drawn. 

The outline of the paper is as follows. In section II we introduce our model and briefiy review the properties 
of the stringy black hole solutions. The solutions fall into four categories: neutral, coloured, magnetically charged 
and electrically charged, depending on the behaviour of the gauge field. For the first three types of solutions, the 



linearized perturbation equations decouple, under an appropriate choice of gauge, into two sectors, corresponding to 
gravitational and sphaleronic perturbations. We first concentrate on the sphaleronic sector, and show in section III 
that there are topological instabilities for both coloured and magnetically charged black holes. We then count the 
number of unstable modes in this sector and find that it equals the number of zeros of the gauge field function of the 
background solution. In the same section, we also consider the stability of the sphaleronic sector of the neutral black 
holes. In section IV, we comment on the gravitational sector for all four types of black holes, appealing to catastrophe 
theory and continuity arguments. The magnetically charged and coloured black holes have instabilities in this sector 
due to the presence of the non-Abelian gauge field, whilst the bulk of the electrically charged and neutral black holes 
are shown to be stable. Section V is devoted to the conclusions derived from our analysis. 



II. STRINGY COLOURED BLACK HOLES 

We consider the following Lagrangian describing stringy corrections to the SU(2) Einstein- Yang-Mills model |12] 

^ = ^ + \^^^^'^+^ ^f^'^GB - F'-^'^n^) , (2.1) 

where cj) is the dilaton field and TZq^ — Rf^^paR^'^'"^ — iR^^R^'^ + R^ the higher-derivative Gauss-Bonnet term. In 
equation ( |2.1[ ) we have set the gravitational coupling k^ — SttG equal to unity, so that a' /g^ is the single coupling 
constant. In effective string theory, which is the situation in which we are interested, the constant f3 is equal to unity, 
and from now on we shall fix f3 to have this value. It is known that, in this case, there are no particle-like solutions 
to the field equations jlj]. In this paper, we are concerned only with the black hole solutions, and shall not consider 
the particle-like solutions which exist for other values of /3. However, much of our analysis applies equally well for 
all values of (3, and in particular to the case /3 = 0, which corresponds to the Einstein- Yang-Mills-Dilaton (EYMD) 
black holes considered in [^[Q]. If the coupling constant is fixed (for numerical computations the value a' /g^ = 1 is 
convenient) then there is only one remaining parameter, namely the event horizon radius r^. 

We consider the general spherically symmetric ansatz for the SU{2) non-Abelian gauge potential ||2(](| 

A = aoTr dt + b-fr dr + [vfo — (1 + w)T^p\ dO + [(1 + id)t0 + vf^p] sinOdip , (2.2) 



with TV = T.Cr, where fi, i = 1, 2, 3 are the usual Pauli matrices. This ansatz for the gauge potential (2.2) does not 
completely fix the gauge. There is still freedom to make unitary transformations of the form 



TAT-^ + TdT-^ (2.3) 



where 



r = exp [k{r, t)Tr] , (2.4) 

under which the gauge potential components transform according to 



/ oo \ / ao - k \ 
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w 



b-k' 

w cos k ^ V sin k 

\ V cos k + w sin k I 



(2.5) 



where here, and in the rest of the paper, we use ' to denote derivative with respect to r, and ' to denote derivative 
with respect to t. We shall make use of this gauge freedom to choose the gauge for which the perturbation equations 
take their simplest form. 

We also consider the following spherically symmetric line-element for the spacetime background 

ds^ = -e^ dt^ -f e^ dr"^ + r^ {d0^ + sin^ edcp"^) . (2.6) 

In the above two ansatzes for the non-Abelian gauge field and the line-element, F, A, oq, &, v and w depend on 
t and r. For the equilibrium solutions, b and v vanish identically, and all quantities depend on the co-ordinate r 
only. Solutions are found by numerical integration, requiring that there is a regular event horizon aX r — rh and that 



the geometry is asymptotically flat, with no singularities outside the event horizon. The requirement of asymptotic 
flatness places the following constraints on the matter fields as r — > oo: the dilaton (p ^ 0, whilst the gauge field 
function ap ^ 0. Finally, w ^ ±1 for the coloured black holes with vanishing charge at infinity, and w ^ for the 
globally magnetically charged solutions. The types of solutions fall into four categories, depending on the form of the 
gauge field: 

1. Neutral black holes, for which w = ±1 and oq = 0. These solutions have a vanishing gauge field strength and 
were discussed in [ p^|Jl7| . 

2. Magnetically charged black holes jWhere oq = and w = 0. These have a fixed magnetic charge of unity once 
the coupling constants are fixed 13|. 



3. Coloured black holes, where ao = and w varies. These are the stringy sphalerons described in |T2|. The 
function w has at least one zero, and the solutions are characterized by the number of nodes of w. 

4. Electrically charged black holes, where oq does not vanish. The equa tions describing these black holes are found 
from the general field equations (derived by varying the action (2.1)) by setting c = oqw = and then w = ±1. 
The form of the function ao (r) is fixed by the field equations to satisfy 



z[, = e(r+^'/2e 



.Q 



(2.7) 



where the electric charge Q of the black hole can be varied, provided a naked singularity is not formed p3| . 

The form of the metric and dilaton fields are qualitatively the same for all four types of solutions, with the dilaton field 
monotonically decreasing to its asymptotic value and the metric functions interpolating between their horizon and 
asymptotic values at infinity. Regular black hole solutions of each type exist for all values of rh above a critical value 
at which point a naked singularity is formed. In the case of dilatonic and coloured black holes with a Gauss-Bonnet 
term the above constraint takes the form M.Ea 



9' 



< 



V6' 



(2.8) 



The first three types of solution are easily studied within the same algebraic framework, simply choosing w as 
appropriate for each case. Since the electrically charged black holes have a non-zero function aq, they will have to be 
considered separately in the stability analysis. 



III. SPHALERONIC INSTABILITIES 



In this section we consider the magnetically charged and coloured black holes and investigate the existence of 
sphaleronic instabilities under small, bounded, spacetime-dependent perturbations. We use the temporal gauge and 
set ao — 0. The advantage of using the temporal gauge is that the system of perturbation equations decouples into two 
sectors, the gravitational sector consisting of ST, SA, S(j) and Sw, and the sphaleronic sector which comprises Sb and 
Sv. For the electrically charged black holes, the most useful choice of gauge is not immediately apparent. This will be 
considered in section IV, where it is shown that the electrically charged black holes effectively have only gravitational 
sector perturbations. Here, we concentrate on the study of the sphaleronic sector of the coloured and magnetically 
charged black holes leaving the stability analysis of their gravitational sector also for section IV. Our work in this 
section will also be applicable to the stability analysis of neutral black holes under small, gauge perturbations, and 
the corresponding conclusions will be drawn simply by setting w = ±1 in the following. The analysis of subsection 
IIIB will then confirm that the neutral black holes do not have any instabilities in this sector. 

For the sphaleronic sector, and the coloured and magnetically charged black holes, we have the following perturba- 
tion equations, where we have considered periodic perturbations of the form 6V{r, t) — 6V{r)e^''*: 



T' A' 



Sb' + Sb{(j)' + -- — - — ]+—:^wS. 



2e/ 



= 0; 



(3.1) 



and 



r-A / Hbt, Hbu \ _ 2 1 r'e-^ \ f Sb 



"'-"^^.-^J'-'io i (Zl («) 



where 

Hbb = w Sb 



Tibi, = w 5v — w 5v 



Uyb = {w 6b)' + w((l)' + ^ ^ j 6b + w' 6b 



u , r'-A'\ .,, , e 



Hu^ = Siy" +U' + — - — 1 ^t'' + ^(1 - w')6iy. (3.3) 

These equations depend on the dilaton field (f>, however, they have no dependence on the Gauss-Bonnet curvature 
term, except through the static solutions. Therefore, our analysis applies equally well to the case of the Einstein- Yang- 
Mills-Dilaton (EYMD) black holes. In the following subsections, we shall first of all prove the existence of unstable 
modes for both the magnetically charged and coloured black holes, and then proceed to count the number of modes 
of instability. Our analysis will not depend on the details of the equilibrium solutions, only on properties such as the 
number of nodes of the gauge field. In this sense the instabilities are "topological" . 



A. Existence of instabilities 

In this subsection we shall employ a variational method to show that there are sphaleronic instabilities for both the 



coloured and magnetically charged black holes. The sphaleronic sector perturbation equations (3.2) take the form 



»it)-''^(t)' M 



The operator Ti. must be self-adjoint with respect to a suitable inner product on the space of functions 'f 



dv 



and the operator A is required to be positive definite, so that (^|yl|^) > for all non-zero ^P, using the same inner 
product. The variational method has been applied successfully to various systems involving non-Abelian gauge fields 
(see, for example, M), and involves defining the following functional: 



2r,T,i _ (*l^l*) _ W 



^ [*] = M4f*y = M ^'-'^ 



for any trial function \1/. The lowest eigenvalue of the system ( |3.2| ) gives a lower bound for this functional. Therefore, 
there are negative eigenvalues a'^ (which correspond to unstable modes) if we can find any function \1/ which satisfies 
(7^[^] < 0, with (A) < oo. The advantage of this approach is that it is easier to find trial functions which satisfy 
these criteria than it is to find eigenfunctions, which often involves numerical analysis. The disadvantage is that we 
do not obtain precise information about the number or magnitude of the negative eigenvalues. In this subsection we 
are interested in showing that the presence of the dilaton and Gauss-Bonnet term in our model is not sufficient to 
render the gauge field hair topologically stable. We shall return to the question of the number of unstable modes in 
the next subsection. 

The inner product of two trial functions is defined as follows: 

/>OC 

(*|$) = / e'^e-^^-^'>/^^<^dr, (3.6) 



and with respect to this inner product the operator Ti. is self-adjoint, whilst A is positive definite. This inner product 
(3.6) is slightly different from the one usually employed (see, for example, 0]) due to the e'^ term. This term is crucial 



if Ti. is to be self-adjoint. We consider the following trial functions: 

6b = -w'Zk{r) 

6iy^iw^-l)Zkir), (3.7) 



which ar e th e same as those used in the Einstein- Yang-Mills-Higgs case [Q (these are appropriate due to our inner 
product (3_^), which effectively absorbs all the (j) dependence in this sector), fn order to define the functions Zk, first 
introduce a new co-ordinate p: 



dp 
dr 



g0g-(r-A)/2_ 



(3.8) 



This is not the usual "tortoise" co-ordinate because of the (j) dependence. Then we define the sequence of functions 
Zfc(p) by ^ 



Zk{p) = z(£); fc=f,2,... 

where Z[p) is an even function which is equal to unity for p G [0, C], vanishes for p > C -f 1, and satisfies 

forpe [C,C+1], 



^ dZ 
-£>< — <0 
dp 



(3.9) 



(3.10) 



with C and D arbitrary positive constants. 

With these trial functions, after a lengthy calculation, we have 



+ 2{w^-lf] 



/•OO 

{n)=- dr2e^e(r-A)/2j7 

'dr2e^e^^-^'>/^J{l-Zl) 



+ dr2eVr-A)/2(^2_ 1)2^^2^ 



where 



J ^w'^ + -^{w^ - 1)2 > 0, 



(3.11) 



(3.12) 



and all boundary terms have vanished due to the definition of the functions Zk- Immediately it can be seen that the 
expectation values of A and Ti. are fi nite for each value of k, and the expectation value of A is positive, as expected. 
The second and third terms in (3.11) converge to zero as k —^ oo, and hence choosing k sufficiently large, we obtain 
a negative expectation value for H. We conclude that the system of perturbation equations (3.2) has at least one 
negative mode, and hence that the black holes are unstable in the sphaleronic sector. 



The other equation in this sector (3.1) is known as the Gauss constr aint. However, we do not need to consider it 
here, since it is automatically satisfied by eigenfunctions of the system (3.2). To see this, consider the "pure gauge" 
functions given by 



Sb = T', 



6v 



-wT 



(3.13) 



where T is an arbitrary (differentiable) function of r. The perturbations ( t3.13|) satisfy (3^) with eigenvalue cr^ = 
for any T . Therefore, if (^6, bv) are eigenfunctions with eigenvalue a^ ^ 0, they must be orthogonal to the "pure 
gauge" modes: 






A 



5b 
5v 



= 0. 



Performing an integration by parts, we obtain 

/•OO 

/ dre'^e-^(r+A)^2_^g^0^ 

•Irh 



(3.14) 



(3.15) 



where Q is the expression in brackets in the Gauss constraint (3.1). Since J- is arbitrary, and the integrand is assumed 
to be continuous, the Gauss constraint must be satisfied by the perturbations 6b and Sv. 

We note that the analysis of this subsection appUes equall y w ell to the coloured black holes of ||l^ as well as the 
magnetically charged black holes which occur when w = [Q. The instability of the magnetically charged black 
holes may be understood by considering them as coloured black holes in the limit in which the nu mbe r of nodes of the 
gauge field function w goes to infinity. Therefore, they correspond to saddle points of the action (2J) due to the fact 
that the gauge group is non-Abelian. In addition, since the sphaleronic sector perturbation equations do not depend 
on the Gauss-Bonnet term, and the proof of instability depends not on the details of the equilibrium solutions, but 
only on their global features (such as the existence of the event horizon) the analysis also appli es to the EYMD black 
holes Q. We emphasize that the calculations necessary to obtain the expectation value ( 3.11 ) made use only of the 
static equilibrium equation for w (which has the same form for both the Gauss-Bonnet and EYMD models) and not 
that for (f), which docs involve the Gauss-Bonnet term. Therefore the EYMD black holes are also unstable in the 
sphaleronic sector. This is what would be anticipated, but previous studies of the stability of the EYMD system have 
concentrated only on the gravitational sector . Therefore the presence of the dilaton (either with or without an 
associated Gauss-Bonnet term) is not sufficient to remove the topological instabilities of the Yang-Mills field. So far, 
the only way to do this is to introduce a negative cosmological constant |1^ . 



B. Counting the number of unstable modes 



Having shown that there exist instabilities in the sphaleronic sector, we now count the number of unstable modes. 
This subsection extends the method of Q , which counts the number of sphaleronic instabilities of Einstein- Yang-Mills 
(EYM) black holes, to systems involving a dilaton field and Gauss-Bonnet curvature term. 

Firstly, we define the usual "tortoise" co-ordinate r* (compare p (pM)- 



dr* 
dr 



-(r-A)/2 



(3.16) 



and introduce the quantities 



2 



^^!le-(r+^)/V5&, 7^ = ^. 



(3.17) 



Then the perturbed equations for the sphaleronic sector (p^) can be written in the more compact form 



dr* 



a x = w ^ X + e"^ [w -r^^v - -r-6v 



dw 



dr* 



dr* 



(3.18) 



2 d) c ^ 

-a e^w ov — — — 
dr* 



dw 



w 7 x + s [ w — — 6v — — — Sv 



dr* 



dr* 



The first of these equations is the Gauss constraint which we assume to hold even in the case that cr = 0. In that case, 
we can easily prove that the third equation is a direct consequence of the first two and it will be ignored hereafter. 
In terms of the new function 



1 

w 



(3.19) 



the first two perturbation equations in (3.18) can be rearranged to give 



dr* 



-f ZYi(r*)u = cr^u, 



(3.20) 



where 



"■("•' = r'(' + -'' + ^(l^ 



2 dw d(b 1 

— -f - , 

w dr* dr* 4 V dr* 



1 d^ 



2dr* 



(3.21) 



This is a standard Schrodinger equation, but the potential Ui is not regular because w has zeros. Therefore we need 
to use the method of ||6| to map this equation to a "dual" Schrodinger equation which will have a regular potential. 



The "pure gauge" modes (3.13) give a solution of (p-2(\) when a — 0, namely 



wy^ dr* 

where, in order that the Gauss constraint is satisfied (even though cr = 0) the arbitrary function T must satisfy the 
differential equation 

d fe'^ dT\ ., o_ 
-—[—-—] ^e'l'w^T. (3.23) 

dr* V7^ dr* J ^ ' 



The function mq can be used to factorize the operator acting on u in (3.2C): 



\ar* uq dr* J \ar* uq dr* J dr*'' uq dr*'' 

By defining ijj — Q^u and applying Q^ again, we obtain the "dual" eigenvalue equation 

Q-Q+i; = -^+U2{r*)i'^a^^, (3.25) 

where 

„,(..) = 1 y (3„^ - 1, + 2 ^(^^Y) + ii!* + 1 f i^V , (3^26) 



2 ' ^ ^ dr*^ '2 dr*^ A\dr 

In the above, F is a function of r* defined as 



. dr 
and satisfying the equation 



Y = -^'^{^] (3-27) 



dr* -^^ + -^^^ + ^^- (3.28) 



The equation (3.25|) is another standard Schrodinger equation, but now we have a potential which is regular everywhere 



outside the event horizon, provided that K is a regular, well-defined, bounded function for all r* . Assuming for the 



moment that this is the case (we shall return to this issue shortly), the solution iJjq of (3.25) when a = will then 
have the same number of zeros as there are eigenfunctions with a^ < 0, that is, the number of unstable modes. The 
zero energy solution iJjq satisfies the equation 

d 1 duo\ , f d 2,A 1 '^^ 1 



dr* uo dr* I \ dr* w dr* 2 dr* 



leading to the solution 



^0 = u; 6-^/2 g^p / w^Ydr*]. (3.30) 



The function iJjq then has the same number of zeros as w if the last factor on the right-hand-side of ( |3.30D is regular for 
all r* . Assuming, in addition, that i/'o is normalizable, this leads us to conclude that the number of unstable modes 
in the sphaleronic sector equals the number of zeros of the gauge field function w for the coloured stringy black holes. 

At this stage we need to check our assumptions, namely that the function Y , satisfying the differential equation 
( |3.28 ), is regular for all r* and vanishes sufficiently quickly in the asymptotic regimes r* -^ ±oo so that the integral 



in ( |3.30 ) is finite everywher e. W e also require that the function i/>o given by (3.30) is a normahzable eigenfunction. 
Firstly, from the definition ( |3.27 ) of Y, it can be seen that Y will be a regular function as lon g as J- is regular, and 



^p- 7^ 0. The behaviour of JF is most easily found by considering the differential equation (3.23) in terms of the radial 
co-ordinate 



J^" 



r 



y 



T' 



2e^w^ 



■T = 0. 



(3.31) 



This equation has regular singular points at r = r^,oo. Near r = rh, t he sta ndard Frobenius method reveals that 
either JF ~ 0(1), or JF ~ 0{r — rh) as r ^ r^. We consider the solution of (3.31]) which has the behaviour JF ~ (r — rh). 



Since ( 3.31 ) has only one other singular point, at infinity, this solution can be extended to a solution regular for all r. 
Then, as r ^ oo, applying the Frobenius method again shows that either T ~ 0{r) or JF ^ 0{r~^). It mu st be the 
case that T ^ r^ since we can show that T cannot vanish at both r = rh and as r ^ cx), as follows. From (3.23) we 
have 



V— f-— 1 dr* 
dr* \ 7^ dr* _ 



Te'P dT 



ry2 ^j.* 



^ e'P /dT 
r* 7^ \dr* 



dr 



(3.32) 



Taking Tq — > — oo (corresponding to r — > r^), r* — > oo (r — + oo), we obtain a contradiction if J^ — > for both limits (if 
J- is not identically zero). Similarly, if dT /dr* ~ for r* — r^, and taking Tq — > — oo, we also obtain a contradiction, 
which means that dT/dr* cannot be ze ro for r* G (—00,00). Therefore y is a regular function of r* for all r*. In 
addition, the differential equation ( ^.2^ ) shows that 



Y 



as r* 



-00, 



Y ^^0 as r* 



(3.33) 



This means that |^oP ^ (^*) ^ ^^ i"* ^^ ±00, and so ipo is a normahzable wave function. Therefore the conditions 
necessary for our conclusion, that there are as many negative modes as zeros of the function w, to hold are satisfied. 

The analysis of this subsection once again reveals the topological nature of the instabilities, since the precise details 
of the equilibrium solutions were not needed, but only general properties such as the behaviour in the asymptotic 
regions and the number of nodes of the gauge function. In particular, the only equilibrium field equation used in the 
calculations is that for w, which docs not depend explicitly on the Gauss-Bonnet term. Therefore, the result holds 
for all values of (3, so both the coloured black holes with a Gauss-Bonnet curvature term and the EYMD solutions are 
covered. The number of unstable modes is exactly the same as for the EYM black holes j^] , so the stringy corrections 
make no difference to the topological instabilities. 

So far in this subsection we have been concerned with the coloured black holes, since we have assumed implicitly 
that w does not vanish iden tically. For the magnetically charged black holes, where w = 0, the sphaleronic sector 
perturbation equations (3.2) reduce to Sb = (implying that the Gauss constraint is satisfied) and 



2. d^ , 



d(j) d 
dr* dr* 



5v + -^5v = 0. 



^'^ + ik-Us{r*)i = 0. 



This can be cast into the form of a standard Schrodinger equation for the variable ^ = e'^/'^6v: 

dr*^ 
where the potential is 

1 d^<h 1 



Us{r* 



2 dr*'^ 4 V dr 



2 r 
e 



(3.34) 



(3.35) 



(3.36) 



In the next section we shall find that a Schrodinger equation with this potential also governs the gravitational sector 
perturbations of t he ma gnetically charged black holes. We shall be able to appeal to catastrophe theory to show 
that this equation (B.35) has an infinite number of unstable modes, so that the magnetically charged black holes have 
infinitely many sphaleronic instabilities. This result was anticipated from regarding the magnetically charged black 
holes as the limit of coloured black holes in which the number of zeros of w goes to infinity. 

As we mentioned in the beginning of this section, our analysis is also applicable to the stability analysis of the 
sphaleronic sector of the neutral black holes. Although the corresponding background solutions have no ga uge field, 
a sphaleronic sector arises when one applies small gauge perturbations to the system. From equation (3.3C), we can 
easily see that, in this case, the function ipo has no nodes since, by definition, w'^ = 1 everywhere. This confirms the 
absence of unstable modes in the sphaleronic sector of the neutral black holes and demonstrates the stability of this 
family of solutions even under gauge perturbations. 



IV. GRAVITATIONAL SECTOR PERTURBATIONS 



We now turn to the gravitational sector perturbations. The perturbation equations in this sector, as shall be seen 
below, are extremely unattractive (in [[l7| a great deal of complex numerical work was necessary) and so we shall 
appeal to the catastrophe theory analysis of M^. In |13| the authors investigated the properties of the static solutions 
with varying horizon radius rh, fixing the values of the coupling constants o! j g^ and /3. In this section we shall firstly 
consider the gravitational sector of the coloured and magnetically charged black holes. Then, we shall confirm that 
the addition of the perturbation bw in the gravitational sector of the neutral black holes does not introduce any 
instabilities. Finally, we will examine the stability of the electrically charged black holes, which need to be studied 
separately because ag 7^ for the equilibrium solutions. 



A. Coloured and magnetically charged black holes 

The gravitational sector perturbation equations for the neutral, coloured and magnetically charged black holes using 
the temporal gauge take the form: 
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The presence of the Gauss-Bonnet term means that these perturbation equations are considerably more complex than 
those of the EYMD system [^. However, the gauge field does not make them particularly more complicated than for 
the neutral black holes with the Gauss-Bonnet term |1^ . 



We can immediately reduce the number of equations by one, by integrating equation (4.1c) with respect to time to 
give 
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(4.2) 



where /z(r ) is a n arbitrary function of r. Differentiating the above equation with respect to r and comparing it with 
equation (4.1c) by using, at the same time, the time- independent equations of motion, a lengthy calculation gives the 
following differential equation for [i(r) 



m'W + mW 



^(r'-A') + i 
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(4.3) 



When integrated, this equation yields [i(r) ex e^ ^^''"^ /r which goes to infinity when r ^ r^ in contradiction to our 
assumption of small, bounded perturbations. As a result, the only acceptable solution for the function ^{r) is the 
trivial one, /^ = 0. 

For both the coloured and magnetically charged black holes, the behaviour of the static solutions as Vh is varied is 
qualitatively the same. Static configurations which solve the field equations exist for all r/j above a certain value, at 
which point a naked singularity forms. This is in accord with the result that there are no particle-like solutions in 
this theory [|l9[ . There is no critical point in the graph of black hole mass against horizon radius jl^ , and therefore 
catastrophe theory tells us that the stability of the solutions does not alter as we vary r^ . This result can be considered 
from another perspective. 

The system of equations (4.1) could be rewritten as two coupled linear equations for the perturbations 5(j) and 5w 
by eliminating 5K and 6T . (This would be a long and tedious calculation, see p^ for the corresponding computation 
in the neutral case.) Assuming that the equilibrium solutions are analytic in the parameter r/j, it would be possible 
to invoke functional analysis theorems to show that the negative eigenvalues a^ of the system must also be analytic in 
rh- Therefore the number of negative eigenvalues (corresponding to the number of unstable modes) can only change 
when an eigenvalue passes through the value zero. A perturbation changes the mass of the black hole by an amount 
proportional to the eigenvalue. Therefore a zero mode does not change the mass of the black hole. Furthermore, a zero 
mode corresponds to a time-independent perturbation, in other words, a small, static perturbation of the equilibrium 
solutions. Such a perturbation exists only if there are two equilibrium solutions arbitrarily close together for the same 
value of the black hole mass, that is, at a critical point. Therefore, when there is no critical point (as is the case 
for the coloured and magnetically charged black holes), the stability does not change as we vary the parameter ru- 
For the neutral black holes, there are two branches of solutions [Q, the upper branch extending to arbitrarily large 
rh- The analysis below would then apply to this upper branch, and we shall below obtain agreement with the known 
result that this branch of solutions is linearly stable jl^ . 

Since there is no upper bound on the value of r^, we shall consider the perturbation equations as r/j — > oo, in which 
case they will take a particularly simple form. From this we shall be able to appeal to the catastrophe theory analysis 
to deduce the stability of the black holes, whatever their horizon radius. First, introduce the following dimensionless 
variables: 
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In that case, all the terms proportional to a' become of 0(l/r^) while all the other terms are of 0(1). However, 
we cannot yet take the limit r/i — > oo since the Gauss-Bonnet terms contain derivatives of the metric function s tha t 
diverge at the limit r — > r^. In order to resolve this, we use, once again, the tortoise coordinate transformation (3.1(;) 
which now connects f with f*. Then, the perturbed equations take the form 
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As rh —^ oo, in accordance with the no-hair theorem, the geometry becomes that of a Schwarzschild black hole with 
(j) = const and the Yang-Mills field superimposed on this background. As a result, both of the perturbations 5A and 
5T vanish while the equation for 5(1) reduces to the following form 



d^X e(r-A)/2 / dT dA 



df*^ f V df* df* 
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where 



A = ^0exp 



=(r-A)/2 



■df* 



(4.7) 



and we have considered periodic perturbations. The above equation takes the form of a Schrodinger-like differential 
equation with a potential which is everywhere regular as well as positive and vanishes as f* — > ±oo. We, thus, 
may conclude that the sub-sector of the dilaton and metric perturbations is still characterized by the absence of any 
unstable modes. Now we turn to the perturbation equation of the Yang-Mills function. By implementing the results 
derived from the above analysis, this takes the form 



-rr^Sw — IT (1 ~ 3w'^) Sw = a'^5w . 
dr*^ r2 ^ ' 



(4.8) 



This equation is exactly the same, to leading order in 1/r^, as the perturbation equation for the gauge field in Einstein- 
Yang-Mills theory without a dilaton or Gauss-Bonnet term. This was to be expected since the EYM equations also 
decouple in the limit r/j — > oo to give a Schwarzschild geometry with a gauge field on this background. It is known 
that the EYM system possesses instabilities in this sector (see, for example, g). The above equation for 6w is in 
agreement with this result. As the background gauge function w oscillates around its zero value p2|, a potential well is 
formed in the region where w'^ < 1/3, thus, leading to the existence of bound states in the corresponding Schrodinger 
equation. We therefore conclude that the stringy coloured black holes are also unstable and their instabilities are, 
once again, associated with the existence of nodes of the background gauge function w. It is worth noting that, in the 
case of coloured black holes arising in the presence of a negative cosmological constant in the theory, solutions with 
zero number of nodes do exist and they were proven to be linearly stable in both perturbation sectors |11[ . 

That the magnetically charged black holes also possess unstable modes in this sector may be deduced, as in the 
sphaleronic sector, from the fact that they are the limit of coloured black holes in which the number of zeros of the 
gauge field goes to infinit y. Fo r the magnetically charged black holes, which follow if we set w = 0, the potential in 
the Schrodinger equation (|4.8|) is everywhere negative and goes to zero as f* -^ zLoo: 
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uin = -^- 



(4.9) 



then shows that there are an infinite 



The standard estimate for the number of bound states — J^ \/—U{f*) df* 
number of unstable modes in this case. 

This analysis for very large horizon radius has confirmed what might have been anticipated, namely that the 
instabilities are due to the presence of the gauge field. In addition, it is known that the EYMD black holes possess 
gravitational instabilities H, which is in agreement with our results, assuming that stability does not change as /3 
varies. This is a reasonable assumption, since the behaviour of the field equations as r/j ^ cx) is the same for all p. 

As explained earlier in this subsection, catastrophe theory tells us that the stability of the black holes does not 
change as we vary r/j. The results above, derived in the limit r/j — > oo, can therefore be extended to arbitrary r/j 
for which black hole solutions exist. In particular, the coloured black holes will be unstable in this sector, and the 
magnetically charged black holes will have an infinite number of unstable gravitational modes. 

We now exploit the result that the magnetically charged black holes have infinitely many unstable modes in the 
gravitational sector, irrespective of the value of the horizon radius, to infer that they also have infinitely many unstable 
modes in the sphaleronic sector. For the magnetically charged black holes, and any value of the horizon radius rh, 
the perturbation equation for Sw (4.1b) decouples from the other equations in the gravitational sector and has the 
simplified form 
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dr*^ 
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where ^ = e'^^^Sw and 



^^ ' 2dr* 
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(4.11) 



This potential is exactly the same as that arising in the sphaleronic sector perturbations of the magnetically charged 
black holes ( |3.36 ). Therefore, since the Schrodinger equation (4.10) has infinitely many negative eigenvalues, the 
magnetically charged black holes also have infinitely many unstable modes in the sphaleronic sector, as asserted in 
the previous section. 

For the neutral black holes, in the limit of large horizon radius, the potential in the perturbation equation for 6w 



( [4. 8]) , after setting w — ±1, reduces to 



U{f* 



2e' 



(4.12) 



Since this potential is everywhere positive and tends to zero in the asymptotic regions, we can conclude that there are 
no unstable modes. Therefore the neutral black holes are stable as we vary r^ (we emphasize that this applies only 
to the upper branch of solutions in |13[ , namely the branch of solutions which extends to arbitrarily large horizon 
radius and which includes the bulk of the background black hole solutions 1 18|). In the next subsection we shall show, 
by another method, that the equation for Sw and general r^ has no unstable modes. This result will be useful in 
studying the stability of the electrically charged black holes, to which we now turn. 



B. Electrically charged black holes 

We now turn to the electrically charged black holes, which need to be considered separately from the other cases 
because aq ^ for the equilibrium solutions. As in the previous subsection, we first consider the linearized perturbation 
equations involving the metric perturbation (5A: 
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Integrating (4.13b) gives, as in the previous subsection, a function ^{r) such that 
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Since the left-hand-side of equation ( 4.15| ) depends only on r and not on t, the same must be true of the combination 
of perturbations in the brackets on the right-hand-side. For periodic perturbations, it follows that 



Sa'o - a^ (50 + ^a'^{6r + SA) = 0, 



(4.16) 



in which case, exactly as for th e coloured and magnetically charged black holes, /i(r) = 0. It should be stressed at 
this point that equation (4.16) does not in fact correspond to a choice of gauge. According to the residual gauge 
freedom of the gauge potential (pj) the quantity Sb — Suq is invariant, so we cannot choose its form. Equation (4.16) 



therefore represents a constraint on the gauge potential perturbations, necessary in order that the equations (4.12) 
are consistent. When ag = (for the neutral, coloured and magnetically charged black holes) the two equations for 
(5A are automatically consistent, regardless of the behaviour of w. The remaining gauge freedom could be used to set 
(5ao = in the electrically charged case, as for the other types of black hole solution. 

For the construction of the Yang-Mills perturbation equations, we first define the quantity c = aow. For the 
equilibrium black holes, c is set to vanish identically, although w = ±1 and ao ^ ||l^. Therefore, we regard c as a 
separate entity, with perturbation Sc, which is independent of 6w and Sao. The Yang-Mills perturbation equations 
then take the following form: 
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where the ± depend on whether w = ±1. With the constraint ( 4.16 ), the Yang-Mills perturbation equations ( 4.17 ) 
simplify greatly, and reduce to the following conditions: 
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and the equation for Sw then also has a simple form 
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(4.19) 



It is rather surprising that this equation has no time dependence. In order to study the space dependence of Sw, we 
are going to consider the above perturbation as the zero eigenfunction of the following eigenvalue problem 
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We note as an aside that ( [l.2C| ) is the equation satisfied by periodic perturbations Sw in the case of neutral equihbrium 
black holes. Here we shall be using ( 4.20| ) as a tool to show that there are no solutions of ( 4.19| ) representing physical 
perturbations. It should be emphasized that ( 4.19 ) does not restrict the time depe nden ce of the perturbation 6w at 
all, although the spatial dependence for each time t will be that of the solution of ( 4.2C| ) when a^ = 0. 

In order to get ( 4.20|) in the form of a Schrodinger equation, it is convenient to define another new co-ordinate, TZ, 
by 
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dr 
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Again, the presence of the e *" term means that this is not the "tortoise" co-ordinate (compare (3^) and ( 3.16 )) 
Then ([4.20D becomes 
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This equation is not quite of the standard Schrodinger form, due to the e^*^ multiplying the a^ . However, since 
e^*^ > 0, the standard theorems still apply (see, for example [p3[), so in particular there can be no negative eigenvalues 
CT^ because the potential 



U 



26^62-^ 



(4.23) 



goes to zero as 7^ — > ±oo (i.e. at the event horizon and at infinity), and C/ > everywhere. Since there are no 
negative eigenvalues for a"^^ this means that the solution of ( J4.22D with cr = (if it exists) can have no zeros. For a 
physical perturbation 5w, we require that 5w -^ Q as TZ -^ ±oo. Since Sw has no zeros, it must be of one sign and 
have at least one maximum (if it is everywhere positive) or at least one minimum (if negative). By ( 1.22| ), (P5w / cTR? 
has the same sign as 5w, since the potential U is positive. Therefore, if 5w > 0, at a stationary point d Sw/dTZ^ > 
also, which means that Sw has a minimum. This is in contradiction with the requirement that 6w tends to z ero in 
the asymptotic regime. A similar argument holds if dw is everywhere negative. The only possible solution of (4.1£) 
is therefore 6w = 0, so that 



Sw = Sv ^ Sb = Sc = 0. 



(4.24) 



We note that the fact that there are no negative eigenvalues a^ of ( [4.20 ) confirms that there are neutral black holes 
with no instabilities, even when embedded in the non-Abelian gauge group. 

For the electrically charged black holes we are left, as with the neutral black holes, with effectively only a gravita- 
tional sector, consisting of 6(j), 6T, i5A, and Sao, the latter being constrained by ( 4.16|) to be given by: 
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We could use the remaining gauge freedom to set da^ = 0, but this will make no difference to our analysis. The 
remaining perturbation equations take the form: 
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where SA is given by (4.14). This could be reduced to a single equation by using (4.14) and (4.25). However, since 



we shall be appealing to catastrophe theory, this does not afford any advantage in the analysis. 



Now that the system has been reduced to a set of coupled equations ( 4.26 ) involving the perturbations of the 
dilaton field and metric functions, we can once again appeal to the catastrophe theory analysis of p3[ . They find two 
branches of electrically charged black hole solutions, one of which extends out to infinite horizon radius. We shall 
focus on this branch of solutions since it is this branch (if any) which can be stable. As in the previous section, we 
consider the perturbation equations in the limit r/i ^ oo, and work to leading order in 1/rf^- Using the same tortoise 
co-ordinate f* as in the previous subsection, we find the equations: 
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(4.27) 



These are e xact ly the same equations (without the equation for Sw) as those obtained for the other types of black 
hole solutions (4.5). In this case, for r/j ^ 1, the geometry reduces to a Reissner-Nordstrom black hole, with = const. 
As in the Schwarzschild case, the perturbations ST and SA vanish identically, and the equation for 50 reduces to (^^) . 
The same argument then shows that there are no instabilities in this sector for r/i 3> 1. The catastrophe theory 
analysis then tells us that the upper branch of electrically charged black holes is stable for all r^ . 

We conclude in this subsection that the upper branch of electrically charged black holes is stable, whereas we have 
already shown that the magnetically charged black holes are infinitely unstable. This may be surprising, especially 
since the gauge field in both cases is essentially Abelian. However, for the magnetically charged black holes, the gauge 
potential is embedded Abelian, in other words, it corresponds to the product of a U(l) gauge potential and a constant 
matrix. This embedding in the non-Abelian gauge group SU(2) gives rise to the infinite number of unstable modes. 
As discussed earlier in this article, the same conclusion can be reached from observing that the magnetically charged 
black holes arise as the limit of coloured black holes in which the number of nodes of the gauge function w goes to 
infinity. This leads to infinitely many modes of instability in the sphaleronic sector. On the other hand, due to the 
construction of the electrically charged black holes (in particular, setting c = gqw = in the field equations although 
neither oq nor w vanish) means that the gauge potential is genuinely C/(l) (without any embedding). Therefore the 
stability of these solutions is not unexpected. In this case there is effectively no sphaleronic sector, so the other 
instabilities present for the coloured and magnetically charged black holes do not arise. 



V. CONCLUSIONS 

In this article, we have considered a generalized, string-inspired theory of gravity that describes the non-minimal 
coupling of a single scalar field, the dilaton, to gravity through the higher-derivative Gauss-Bonnet term as well as to a 
non-Abelian SU(2) gauge field. This theory has been shown in previous work ||T2] to admit regular and asymptotically 
fiat black hole solutions that are characterized by the presence of a non-trivial dilaton and gauge field on the region 
outside the horizon in contradiction with the "no-hair" theorem of the Theory of General Relativity. Nevertheless, the 
hair of these black hole solutions is merely "secondary" in the sense that no new charges can be associated with the 
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aforementioned, non- vanishing fields. The dilatonic black holes that arise in the same framework but in the absence 
of the gauge field |I^ have been already proven |lj to be linearly stable under small, bounded, spacetime-dependent 
perturbations as they correspond to the stable branch of the family of neutral black hole solutions. Then, the question 
of the behaviour of the coloured black holes under the same type of perturbations, and in the presence of the same 
stringy corrections, naturally arises. 

By making an appropriate choice of gauge, the linearized perturbation equations, for the coloured black holes, were 
decoupled into two sectors, the sphaleronic and the gravitational one. In the first sector, the perturbation equations 
resembled those of the Einstein- Yang-Mills-Dilaton theory with the explicit dependence on the Gauss-Bonnet term 
having been eliminated. For the needs of our analysis, well-known methods [^, 7J, that were previously used for the 
stability analysis of black hole solutions arising in the framework of the Einstein- Yang-Mills theory, were extended in 
order to accommodate the dilaton field. Then, the existence of topological instabilities was analytically demonstrated 
by making use of the method of trial functions. The number of unstable modes in this sector was determined by 
mapping the irregular Schrodinger equation of the gauge perturbations to a "dual" regular one and it was found to be 
the same as the number of zeros of the background gauge field. In the gravitational sector, due to the complexity of 
the perturbation equations, continuity arguments based on the results derived from the catastrophe theory analysis 
p3[ allowed us to work in the limit of infinitely large horizon value. Although the sub-system of the metric and 
dilaton perturbations was found to be stable, instabilities attributed once again to the oscillating behaviour of the 
background gauge function around zero were proven to exist. As a result, we may conclude that the accommodation 
of stringy corrections to non-Abelian black holes fails to render these solutions linearly stable. 

The perturbation equations for the magnetically charged and neutral black holes, under the same type of perturba- 
tions, may easily follow from the corresponding ones for the coloured black holes by simply choosing an appropriate 
value for the background gauge function w. This allows us to draw conclusions concerning the stability of these 
two families of solutions in the same framework and by using the same methods as above. Then, the study of the 
sphaleronic sector reveals the existence of an infinite number of unstable modes, for the magnetically charged black 
holes, and the same holds for the gravitational sector. The above result is in accordance with the interpretation of 
this family of solutions as the limit of coloured black holes in which the number of zeros of w goes to infinity. On 
the other hand, no unstable modes are found for the upper branch of the neutral black holes in both sectors which 
confirms the stable character of these solutions not only under metric and scalar perturbations |17|] but even under 
gauge-dependent perturbations. 

Finally, the stability analysis of the electrically charged solutions was conducted although in a different framework 
of perturbation equations due to the different ansatz for the background gauge field. Nevertheless, we were able to 
show that, in order for our perturbation equations to be consistent, a constraint that involves a combination of gauge, 
metric and scalar perturbations must be satisfied. In that case, we have shown that the gauge perturbations are 
completely decoupled from the scalar and metric ones, a feature which facilitates the study of each subsector. The 
gauge perturbations were found to reduce to a single equation for Sw which, however, was shown not to accept any 
solutions representing physical perturbations. On the other hand, by using the same continuity arguments and working 
again in the limit of infinitely large horizon, we proved that no unstable modes arise in the remaining subsector of 
metric and scalar perturbations for the upper branch of the electrically charged black hole solutions. In conclusion, 
the bulk of the electrically charged black holes, in common with the neutral dilatonic black holes, are stable under 
small, spacetime-dependent perturbations. The stability of these two families of solutions can be justified by their 
interpretation as the generalization of the Reissner-Nordstrom and Schwarzschild black holes, respectively, in the 
framework of string theory. 
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